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BEHAVIOR  OF  POROUS  BERYLLIUM  UNDER 
THERMOMECHANICAL  LOADING: 

PART  4. 

CONSTITUTIVE  MODEL  FOR  WAVE  PROPAGATION 


Abstract 


A new  constitutive  model  for  wave 
propagation  in  porous  materials  is  con- 
structed here.  Based  on  similar  models 
proposed  earlier  by  Herrmann  and  Holt, 
the  new  model  features  (1)  inclusion  of 
deviatoric  stresses,  (2)  an  option  for 
porosity-dependent  relaxation  time  for 
pore  closure,  (3)  elastic  reopening  of  the 
pores,  and  (4)  an  improved  plastic  com- 
paction function. 

This  model  has  been  applied  to  plasma- 
sprayed  porous  beryllium  in  the  as- 
sprayed  condition  (~*14%  porosity)  and 

. ju  " r\JL 

after  sintering  (~  10%  porosity).  We  dis- 
cuss techniques  used  to  correlate  data 
from  the  static  and  dynamic  experiments 


which  were  conducted  to  provide  the  sup- 
porting inputs  to  the  modelj"  Wave  pro- 
files generated  from  plate  impact  experi- 
ments were  used  to  test  the  model  under 
different  experimental  conditions.  The 
resulting  good  agreement  between  the 
model  calculations  and  wave  profile  ex- 
periments suggests  that  the  mechanical 
response  of  porous  beryllium  can  be  satis- 
factorily described  within  the  necessary 
accuracy  by  this  model.  We  also  discuss 
two  other  relevant  topics,  namely,  the 
effect  of  an  in-material  gage  on  observed 
wave  profiles,  and  a transient  phenomenon 
affecting  the  wave  profile  at  the  impact 
surface,  j 


1.  Introduction 


When  a thin  plate  impacts  at  high  veloc- 
ity against  a metal  target,  it  generates  a 
stress  pulse  which  travels  through  the 
target  medium.  Response  to  the  stress 
wave  is  slightly  more  complicated  when 
the  medium  is  porous.  First  we  note  that 
the  compressive  side  of  the  pulse  travels 
with  much  lower  velocity  in  the  porous 
medium  than  in  the  corresponding  solid 
medium.  In  the  case  of  10%  porous 
plasma-sprayed  sintered  beryllium,  for 
example,  a 2-GPa  (=  20-kbar)  stress  wave 


attains  a speed  of  0.  35  cm/^s,  which 

represents  only  42%,  of  the  shock  speed  of 

solid  beryllium.  This  dramatic  decrease 

in  shock  wave  velocity  is  a common 

occurrence  in  most  porous  metals  and  is 

associated  with  the  extra  time  taken  by 

the  shock  wave  either  in  closing  the  pores 

or  in  being  scattered  from  them.  The 

second  notable  feature  occurs  on  the 

release  side  of  the  pulse.  Experiments 

on  the  release  wave  velocities  in  porous 
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aluminum  and  beryllium  have  shown 
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that  these  velocities  depend  less  sensi- 
tively on  initial  porosity  and  lie  closer  to 
release  wave  velocities  in  corresponding 
nonporous  samples.  A 2-GPu  stress 
wave  in  1 0"'.  porous  sintered  beryllium 
quoted  in  the  above  example  has  an  initial 
release  wave  speed  of  1.3  cm  ps.  This 
is  nearly  the  same  as  the  release  wave 
speed  of  solid  beryllium.  Both  of  these 
features  imply  that  reopening  of  the  pores 
does  not  take  place  to  a large  extent  and 
that  the  pore  closure  achieved  during 
compaction  is  primarily  accomplished  by 
a plastic  flow  process. 

The  large  difference  between  the  com- 
pressive and  release  wave  velocities 
shows  up  more  dramatically  in  the  atten- 
uation of  a stress  pulse.  The  pulse 
attenuates  very  rapidly  after  propagating 
a relatively  short  distance  through  the 
porous  medium.  This  immediately  sug- 
gests a useful  application.  Namely,  a 
thin  layer  of  porous  material,  coated  over 
a solid  substructure,  can  be  used  to  pro- 
tect the  substructure  against  high-intensity, 
short-duration  shock  waves  arising  from 
either  mechanical  impact  or  other  types 
of  rapid  energy  deposition  (e.g. , sudden 
heating). 

For  such  applications,  it  is  essential 
to  carry  out  computer  simulations  using 
a physical  model  that  can  accurately  pre- 
dict the  dynamic  behavior  of  porous  ma- 
terials. The  present  work  is  addressed 
to  this  problem. 

As  a representative  material  for  the 
porous  metals  we  have  chosen  porous 
beryllium,  whose  dynamic  response  is 
porosity-dependent  up  to  stress  levels 
as  high  as  4 GPa.  Beryllium  remains 
porous  over  a relatively  large  range  of 
pressure,  which  makes  it  useful  for 


testing  tiie  model  at  different  stress 
levels.  I wo  plasma- sprayed  varieties 
were  considered  in  our  modeling  effort  — 
an  "as-sprayed"  beryllium  with  a nomi- 
nal density  of  8(1"; , and  a "sintered" 
beryllium  which  was  densified  to  9CP’<  of 
solid  density  by  a sintering  process, 
experimental  evidence  described  later 
indicates  that  the  two  porous  beryllium 
specimens  show  different  compaction 
and  release  behaviors.  These  features 
were  incorporated  in  our  model. 

A combined  theoretical  and  experi- 
mental study  was  performed  in  which 
structural  characterization  of  the  porous 
beryllium  by  photomicrography,  static 
stress-strain  testing,  and  ultrasonic  data 
provided  information  essential  for  formu- 
lating and  calibrating  the  model,  while 
plate  impact  tests  were  conducted  to  check 
the  predictive  capabilities  of  the  model. 
These  experiments  have  been  carried  out 
at  this  Laboratory,  and  some  of  them 

have  been  reported  separately  and  are  the 
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subjects  of  other  reports  in  this  series. 

Our  model  is  based  on  the  P-a-r  model 

g 

of  Holt  et  al.,  which  is  a rate-dependent 
modification  of  the  P -a  model  developed 
by  Herrmann.^  The  porosity  parameter 
a,  which  is  the  ratio  of  the  porous  to 
solid  volumes,  and  the  relaxation  time  r, 
which  is  the  time  necessary  for  the  pores 
to  close  to  a final  equilibrium  size,  play 
important  roles  in  these  models.  These 
models  develop  a simple  constitutive 
relation  between  the  pressure  and  a.  Thus 
the  difficult  task  of  developing  a compli- 
cated relation  connecting  the  pressure, 
volume,  and  energy  of  a porous  material 
is  bypassed. 

The  difference  in  the  compressive  and 
release  wave  velocities  noted  earlier  is 
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attributed  in  our  model  to  the  change  in 
the  porosity  parameter  o,  while  t can 
account  for  some  persuasive  evidence 
(e.g.,  nonzero  shock  risetime,  disper- 
sion and  attenuation  of  a precursor,  and 
the  possible  presence  of  transient  effects 
in  the  shock  front  when  measurements 
are  taken  at  the  impact  surface)  that 
such  rate-dependent  effects  are  indeed 
present  in  experimentally  measured  wave 
profiles.  Holt's  and  Herrmann's  models 
are,  however,  hydrodynamic  models. 
Therefore,  they  do  not  explain  in  a self- 
consistent  manner  any  effect  (such  as  an 
elastic-plastic  nature  of  release  wave 
profiles)  originating  from  deviatoric 
stresses. 

We  added  the  deviatoric  corrections  to 
these  models.  This  was  done  by  using 
the  porous  bulk-  and  shear-modulus  ex- 
pressions derived  from  slightly  general- 
izing Mackenzie's  expressions1  which 
are  applicable  to  a linear  elastic  medium 
containing  spherical  pores.  These  ex- 
pressions (except  for  a slight  modification) 
are  the  same  as  those  proposed  earlier 
by  Seaman.11  The  present  model  also 
includes  an  improved  static  compaction 
function  that  was  deduced  from  the  hydro- 
static measurements. 

Sintered  and  as-sprayed  beryllium 
behave  differently,  especially  along  initial 
loading  portions.  While  the  extended 
forms  of  Mackenzie's  expressions  are 


suitable  to  describe  the  loading  behavior 
of  sintered  porous  beryllium,  the  same 
procedure  cannot  be  used  to  describe  the 
loading  data  of  the  as- sprayed  porous 
beryllium.  This  is  due  to  a relatively 
large  number  of  nonspherical  pores  and 
irregularly  shaped  cracks  which  separate 
different  beryllium  grains  in  the  as- 
sprayed  material.  For  as-sprayed  porous 
beryllium,  therefore,  we  assume  that  the 
loading  occurs  plastically  from  the  begin- 
ning of  compaction.  For  this  purpose, 
we  introduce  a porosity-dependent  relaxa- 
tion time  t whose  expression  was  derived 
from  a spherical  model  calculation  for  an 
elastic-viscoplastic  material  by 


Cristescu1"  and  Holt  et  al.1^* 


The  use 


of  such  an  expression  gives  a wave  front 
which  is  sharply  rising  at  the  foot  but 
later  becomes  dispersive  at  the  peak,  a 
b^'-avior  much  in  evidence  in  the  wave 
p routes. 

1 4 

A computer  program  describing  the 

model  was  included  as  a constitutive 

model  in  the  one-dimensional  hydrody- 
1 5 

namic  code  KO  and  was  applied  to  plate 
impact  experiments,  where  accurate  data 
on  full  and  attenuated  wave  profiles  were 
available.  This  comparison  provides  one 
of  the  severest  tests  of  the  model's  ade- 
quacy. The  satisfactory  results  described 
below  lead  us  to  believe  that  the  mechan- 
ical response  of  porous  beryllium  can  be 
described  adequately  by  our  model. 


2.  Characterization  of  Porous  Beryllium 


Our  porous  model  requires  information  tered  Be  plates  contained  about  1.5  wt% 

from  several  types  of  static  and  dynamic  of  BeO.  The  other  impurities  (Fe,  Al,  C, 

experiments  on  porous  beryllium.  These  Mg,  W)  totaled  0.35%  and  0.20%  each  for 

experiments  were  done  on  porous  beryl-  the  as-sprayed  and  the  sintered  specimens, 

lium  specimens  from  the  same  batches  A mure  complete  description  of  the  two 

in  all  the  test  series  to  avoid  the  common  materials,  including  spectrographic 
problem  in  experiments  on  porous  mate-  analyses,  is  given  in  Ref.  3.  In  the 

rials  of  batch-to-batch  variations,  thus  computations  presented  later,  we  use 

assuring  a consistent  comparison  between  initial  densities  of  1.660  g/cm3  and 

3 

experiments  and  computations.  Several  1.591  g/cm  for  the  sintered  and  the  as- 

groups  at  LLL  contributed  to  this  work:  sprayed  Be  specimens,  respectively, 

photomicrography  and  material  charac-  These  values  agree,  within  the  experi- 

terization  data  was  provided  by  Hanafee  mental  uncertainties,  with  measurements. 

4 

and  Snell,  quasi-static  deformation  data  Figure  1 shows  photomicrographs  of 

by  Schock,  Abey,  and  Duba,  ultrasonic  the  sintered  and  as-sprayed  samples. 

g 

velocity  measurements  by  Green,  and  Comparison  of  the  sintered  and  as- sprayed 

Hugoniots  and  wave  profiles  by  Isbell  and  specimens  shows  that  the  pores  in  the 
3 7 

Horning'  and  Gust.  Some  of  these  data  as-sprayed  specimen  are  more  irregular 


had  to  be  reduced  into  forms  applicable  in  shape,  larger  in  number,  and  have  a 

for  model-building  purposes.  We  give  wider  size  distribution.  There  are  also 

below  an  overview  of  the  experimental  cracklike  structures  in  the  as-sprayed 

results  and  discuss  techniques  used  to  material,  formed  by  colonies  of  larger 

correlate  those  features  that  are  relevant  pores.  Sintering,  in  addition  to  densifying 

to  the  present  work.  the  material  by  about  4%,  reduces  the 

number  of  pores  and  makes  them  less 

2.1  STRUCTURAL  CHARACTERISTICS  irregular  in  shape  and  more  uniform  in 

size.  This  difference  is  taken  into  ac- 
Beryllium  powders  supplied  by  Kawecki  count  in  our  model. 

Berylco  Industries  were  plasma- sprayed  In  all  test  specimens,  the  outer  0.015- 

in  a dry  argon  atmosphere  onto  a rotating  cm  layer  from  the  top  and  bottom  sides  of 
aluminum  mandrel,  which  was  removed  both  sintered  and  as- sprayed  machined 

after  the  plasma  spray.  Thickness  of  the  specimens  was  removed  by  chemical 
finished  plate  ranged  from  0.48  to  0.55  cm  etching.  The  densification  of  surface 

across  the  radius.  One  specimen  (desig-  layers  arises  from  the  machining  process 

nated  as  "sintered")  was  heat-treated  and  can  lead  to  erroneous  experimental 

2 hours  at  117  5°C,  while  the  other  speci-  results  if  not  removed.  Considerable 

men  (designated  as  "as-sprayed"  and  with  twinning  of  the  beryllium  grains  was 

a slightly  different  initial  powder)  had  no  noted  in  the  densified  layers,  the  ef- 

such  heat  treatment.  Chemical  analysis  feet  being  most  predominant  in  the 
showed  that  the  as-sprayed  and  the  sin-  sintered  material. 
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Fig.  1.  Pnotomicrographs  of  plasma-sprayed  porous  beryllium:  (a)  sintered  specimen, 
baked  at  450°C,  (b)  as-sprayed  specimens,  baked  at  450°C.  These  cross  sec- 
tions are  normal  to  the  direction  of  plasma  spray.  (Photographs  taken  by 
Hanafee  and  Snell,  Ref.  4.) 


2.2  STATIC  DEFORMATION 
CHARACTERISTICS 

Figure  2(a)  shows  the  hydrostatic 
pressure- volume  relation  for  the  sintered 
beryllium.  The  portion  of  the  hydrostat 
from  0 to  0.8  GPa  was  obtained  under  true 
hydrostatic  conditions,  in  which  a speci- 
men was  coated  with  semiflexible  epoxy 
and  fitted  with  foil  strain  gages  to  measure 
circumferential  and  axial  strains.  It  was 
then  placed  in  a pressure-transmitting 
fluid  (oil),  and  the  pressure-volume  re- 
lationship was  measured.  Above  0.8  GPa 
the  measurements  were  obtained  by  a 
"quasi-hydrostatic"  method,  in  which 
the  sample  was  surrounded  by  a pressure- 
transmitting  medium  (tin)  and  forces  were 
applied  to  top  and  bottom  by  pistons.  Tin, 
with  its  low  yield  strength,  acts  like  a 
fluid  for  pressures  above  its  yield  stress 
(<0.01  GPa).  Note  that  the  hydrostat 
exhibits  an  elastic  yield  at  approximately 
0.4  GPa.  The  slope  of  the  porous  hydro- 
stat at  the  relatively  high  pressure  of 


4 GPa  is  only  slightly  less  than  that  of 
1 7 

the  solid  hydrostat,  but  a relatively 
large  porosity  (2.2%)  still  remains  at  this 
pressure.  The  residual  porosity  of  about 
2%,  which  appears  to  remain  even  under 
extremely  compressed  states/  is  taken 
into  account  in  our  model. 

The  hydrostat  for  as-sprayed  beryllium 
given  in  Fig.  1 of  Ref.  5 was  obtained  by 
using  the  quasi-hydrostatic  data  above 
1.2  GPa  and  the  hydrostatic  data  below 
0.5  GPa.  In  the  region  between  1.2  and 
0.5  GPa,  which  corresponds  to  a relatively 
large  range  in  the  volume  strain,  "inter- 
polated" values  were  used  to  smoothly 
connect  the  lower  and  upper  portions  of 
the  hydrostat.  At  pressures  above 
2.5  GPa  the  hydrostat  for  as-sprayed 
beryllium  lies  below  that  for  sintered 
beryllium  with  a nearly  constant  volume 
offset  of  about  0.008  cm  /g.  Possible 
errors  in  the  experimental  system  for 
quasi-static  compaction  include  fric- 
tional forces  that  exist  at  interfaces 
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Fig.  2.  Hvdrostats  and  Hugoniots  of  porous  and  solid  beryllium:  (a)  sintered  and 

(b)  as-sprayed  specimens.  Solid  dots  and  triangles  are  experimental  flugoniot 
points  of  sintered  and  as-sprayed  specimens,  respectively.  Hydrostat  for 
sintered  specimen  is  that  of  Schock  et  al.,5, 16  while  hydrostat  for  as-sprayed 
specimen  is  obtained  by  an  interpolation  scheme  described  in  the  text. 


of  the  sample,  the  die,  and  the  piston, 
and  problems  arising  from  the  exper- 
imental necessity  of  stacking  several 
layers  of  porous  beryllium  disks  to 
make  a specimen  thick  enough  to 
test.  Such  errors  are  of  a magnitude 
to  account  for  the  observed  volume 
differences. 

Therefore,  for  the  calculations,  we 
used  a slightly  different  hydrostat 
(F'ig.  2(b))  for  as-sprayed  beryllium,  but 
one  which  is  probably  equally  reliable. 
This  was  drawn  using  the  following  physi- 
cal observations.  First,  we  noted  that 
the  Hugoniots  of  both  the  sintered  and  the 
as-sprayed  specimens  shown  in  Fig.  2(b) 
are  nearly  equal  to  each  other  for 
stresses  larger  than  about  1 GPa.  This 
suggested  to  us  that  variations  in  the 
microstructure  of  beryllium  are  low- 
strength  phenomena,  and  that  they  do  not 


affect  the  bulk  properties  at  higher 
stresses.  We  concluded,  therefore,  that 
the  hydrostats  of  both  as-sprayed  and 
sintered  specimens  would  also  be  close 
to  each  other  if  the  pressure  exceeded  a 
certain  level,  probably  about  1 GPa. 
Additional  insight  was  obtained  from  the 
ultrasonic  measurements.  The  experiment 
which  was  conducted  to  detect  a transmitted 
shear  wave  failed  to  do  so  because  of  the 
shear  wave's  small  magnitude  in  the  case 
of  the  as-sprayed  specimen.  Assuming 
that  the  specimen  is  isotropic,  this  im- 
plies that  both  the  hydrostatic  P-V  data 
and  the  one-dimensional  (1-D)  strain  data 
would  likely  have  nearly  equal  initial 
slopes  and  hence  come  very  close  to  each 
other  at  low  stress  levels.  Our  porous 
hydrostat  in  Fig.  2(b)  is  drawn  to  agree 
with  these  observations.  That  is,  at 
pressures  beyond  1.2  GPa  it  is  forced  to 
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agree  w ith  the  porous  hv  irostut  of  the 
sintered  specimen,  and  uelow  0.2  O Pa  it 
is  forced  to  igree  with  ttie  1-1)  strain 
data.  \\  e should  note,  however,  that 
porous  plasma  sprayed  beryllium  is  not 
perfectly  isotropic  at  a low  stress  level. 

The  hydrostat  and  the  mean  stress  ob- 
tained from  tiie  1-1)  strain  tests  do  not 
coincide  >t  i given  volume  strain,  although 
the  average  of  the  two  agrees  with  tne 
axial  Stress  old  lined  from  1 - D strain 
tests  at  0.2  (IPa  or  less.  At  present,  it 
is  not  clear  how  much  the  anisotropic 
correction  would  affect  the  wave  profiles 
at  low  stress  levels.  We  assume  it  neg- 
ligible in  our  treatment.  As  will  be  seen 
later,  wave  profiles  that  were  obtained 
under  this  assumption  agree  satisfactorily 
with  experimental  wave  profiles. 

Another  static  deformation  experiment 
yielding  information  that  we  incorporated 
in  the  model  was  the  cyclic  unloading  and 
reloading  experiment  under  1-D  strain 
conditions.  When  porous  samples  were 
subjected  to  a cyclic  loading-unloading- 
reloading stress  with  ever-increasing 
load  after  every  cycle,  the  reloading 
portion  of  the  stress-strain  path  returned 
very  close  to  the  original  stress  and 
strain  levels  where  the  unloading  started. 
This  feature  was  incorporated  in  our 
yield-strength  expression. 

2.3  ULTRASONIC  VELOCITIES 

The  values  of  the  sound  velocities  given 
below  were  taken  in  the  direction  perpen- 
dicular to  the  surface  of  the  porous  speci- 
mens. This  is  the  same  direction  in  which 
the  plasma  was  sprayed  in  making  the 
specimen  material,  and  it  is  also  the 
direction  of  stress  wave  propagation  in 
our  plate  impact  experiments. 


Ultrasonic 

velocity  (cm/|is) 

Longitudinal  Shear 

Sintered  Be 

1.15  0.81 

As-sprayed  He 

0.34  -0a 

aToo  small  to  be  measured. 


In  the  case  of  the  sintered  specimen, 
the  longitudinal  velocity  was  taken  at 
10  MHz  and  the  shear  velocity  at  5 MHz. 
The  values  quoted  above  are  the  average 
values  taken  for  two  different  thicknesses. 
Probable  errors  in  these  measurements 
are  ±2%.  For  the  as-sprayed  specimens 
a test  of  the  frequency  dependence  of  the 
longitudinal  velocity  was  made  at  5 MHz 
and  15  MHz,  using  two  plates  cut  from 
the  same  specimen.  The  measured  veloc- 
ities at  these  frequencies  were  0.343  and 
0.338  cm  qs  for  the  first  plate,  and  0.358 
and  0.355  cm  ps  for  the  second  plate. 

This  indicates  no  significant  frequency 
dependence  on  the  sound  speed  over  the 
frequency  range  tested.  Since  the  data 
on  the  two  plates  agree  closely,  we  chose 
the  values  obtained  from  the  first  plate 
in  our  model.  As  we  mentioned  earlier, 
measurement  of  the  shear-wave  velocity 
was  not  successful.  Its  small  magnitude 
was  difficult  to  detect  because  of  mode 
conversion,  dispersion,  and  or  attenuation 
that  resulted  from  interaction  between  the 
shear  wave  and  the  inhomogeneities  in  the 
porous  specimen. 

2.4  SHOCK  WAVE 

CHARACTERISTICS 

Three  series  of  plate  impact  tests 
were  conducted.  The  first  series  was 
designed  to  obtain  Hugoniot  data  over  the 
stress  range  0 to  4 GPa,  while  the  second 
series  was  used  to  measure  profiles  of 


I 


attenuated  and  unattenuated  shock  waves 
at  different  specimen  thicknesses  and  dif- 
ferent initial  stresses.  The  third  series 
of  shock  wave  tests  was  conducted  using 

high  explosives  to  measure  Ilugoniots  to 
3 7 

33.5  GPa."’  Hugoniot  points  were  needed 


and  fused  quartz)  with  known  Hugomots 
had  been  mounted.  With  this  arrangement, 
the  buffers  not  only  reduced  the  noise  on  the 
gage  signal  coming  from  the  abrasion  of  the 
front  surface,  but  also,  used  as  impedance- 
matching materials,  allowed  four  different 


to  estimate  the  deviator  contribution  in 
the  model,  while  shock  wave  profiles  were 
used  to  check  the  adequacy  of  our  model. 

Figure  3(a)  shows  the  impactor-target 
geometry  used  to  measure  Hugoniot  points. 
Four  porous  specimens  were  mounted  on 
the  front  of  a projectile  and  impacted 
directly  onto  x-cut  quartz  piezoelectric 
gage  assemblies  on  which  four  different 
buffers  (PMMA,  tantalum,  x-cut  quartz. 


stresses  to  be  obtained  on  a single  shot. 

In  Table  1,  the  measured  Hugoniot 
stress  (aR)  and  particle  velocity  (U^) 
data  as  well  as  the  values  of  the  volume 
and  shock  velocity  (Ug)  computed  by  using 
crR  and  Up  in  Hugoniot  relations  are  sum- 
marized for  both  the  sintered  and  as- 
sprayed  materials.  In  the  calculation  of 
U and  V for  the  sintered  material,  we 
used  the  value  of  0.2  GPa  for  the  stress 


Table  1.  Experimental  data  and  analytically  fitted  Hugoniot  points  for  sintered  and  as- 
sprayed  porous  beryllium.  Up  = particle  velocity,  Us  = shock  velocity, 

CTp j = Hugoniot  stress,  and  V = specific  volume. 


Experimental  data  for 

sintered  porous  Be 

Analytic  fit. 

Eq.  (la) 

U 

P 

(cm/ps) 

U 

s 

(cm/ps) 

cth 

(GPa) 

V 

(cm  /g) 

cth 

(GPa) 

V 

(cm3/g) 

0.00266 

0.5071a 

0.338 

0.6062a 

0.333 

0.6059a 

0.0071  1 

0.4561 

0.678 

0.6004 

0.653 

0.5997 

0.01107 

0.4372 

0.922 

0.5942 

0.918 

0.5939 

0.01684 

0.4186 

1.234 

0.5840 

1.283 

0.5848 

0.03098 

0.4017 

2.216 

0.5636 

2.170 

0.5623 

0.04435 

0.4190 

3.144 

0.5444 

3.173 

0.5448 

0.05601 

0.4598 

4.350 

0.5353 

4.342 

0.5349 

Experimental  data  for  as-sprayed  porous 

Be 

Analytic  fit. 

Eq.  (lb) 

U 

P 

(cm  ps) 

U 

s 

(cm/ps) 

aU 

(GPa) 

V 

(cm3/  g) 

cth 

(GPa) 

V 

, 3 . 

(cm  , g) 

0.00601 

0.3050a 

0.277 

0.6171a 

0.291 

0.6 177a 

0.01005 

0.3161 

0.525 

0.6109 

0.504 

0.6101 

0.0132 

0.3230 

0.680 

0.6045 

0.677 

0.6044 

0.0212 

0.3568 

1.11 

0.5896 

1.129 

0.5903 

0.0366 

0.3518 

2.06 

0.5651 

2.043 

0.5645 

0.0466 

0.3667 

2.69 

0.5501 

2.700 

0.5504 

0.0573 

0.3978 

3.62 

0.5394 

3.618 

0.5394 

The  Hugoniot  relations  are  used  to  obtain  these  data.  For  sintered  porous  Be,  we 
used  a Hugoniot  elastic  limit  cthel  = 0.2  GPa  and  UREI  = 1.15  cm/ps. 
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IMPACTOR 


TARGET 


QuOrtz 
^ages  , 


Known 
impoc  fors 


Porous  Be 
spec  mens 


Quartz 

Duffers 


Quartz 

gages 


Multi -layer  multiple  carDon-Kapton 
immaterial  goges 'sensitive  element 
is  darkened  area) 


Known 
impoc  tor 


Porous  Be 
specimens 


Impactor-target  geometries  used 
to  obtain  (a)  Hugoniot  points, 

(b)  transmitted  wave  profiles 
taken  with  quartz  gages,  and 

(c)  wave  profiles  taken  with  in- 
material gages. 


(^hfi  I a*-  the  Hugoniot  elastic  limit  with 
a velocity  of  1.15  cm/'/us  determined  from 
the  ultrasonic  measurement.  The  cr^-vs- 
U data  were  fitted  with  the  polynomials 
in  U (cm  ps)  given  below: 

n ° 


Sintered  Be 


(GPa)  = 0.1284  + 79.21  - 8741 


\s-sprayed  Hi- 


a,,  (GPa)  44.891  • 682.41  2 

11  P p 

- 1.782  X 104('  2 
P 

+ 1 .999  X 103C  4 
P 

°H  - 4 GPa.  (lb) 

As  shown  in  Fig.  4(a)  and  Table  1,  these 
expressions  give  accurate  representations 
of  the  experimental  points.  The  deviations 
between  the  experimental  points  and  the 
values  calculated  from  Eq.  (1)  lie  within 
the  experimental  uncertainties  (±3%  in 
stress  and  ±2°1  in  particle  velocity). 
Equation  (1)  was  used  to  calculate  vs  V 
(Fig.  2)  and  Ug  vs  U (Fig.  4(b)).  These 
calculations  are  not  affected  to  a significant 
degree  by  the  choice  of  0-2  GPa. 

The  HEL  was  measured  as  a function  of 
distance  of  propagation  and  was  seen  to 
decrease  from  about  0.4  GPa  at  0.1  cm  to 
about  0.2-0. 3 GPa  for  distances  greater 
than  0.3  cm.  Figure  4(b)  also  shows 
theoretical  longitudinal  sound  speeds  that 
were  calculated  using  Eq.  (lb)  and  the 
bulk  moduli  of  the  porous  beryllium  speci- 
mens given  in  Sec.  4. 

Among  several  impactor-target  geom- 
etries used  in  the  second  plate-impact 
test  series,  we  shall  concern  ourselves 
with  only  two  types  here.  In  the  first 
geometry,  the  gage  with  its  x-cut  quartz 
buffer  was  placed  on  the  back  surface  of 
the  beryllium  target  and  the  transmitted 
compressive  wave  profile  was  measured 
(Fig.  3(b)).  In  the  second  geometry,  the 


The  use  of  the  in- material  gages  is 
advantageous,  since  they  measure  in  situ 
stress  profiles  without  the  distortion 
caused  by  the  major  discontinuity  pre- 
sented by  a "rear  surface"  gage  or  a 
free  rear  surface.  Such  techniques  can 
reflect  large  compressive  or  release 
waves  which  affect  the  remainder  of  the 
wave  in  the  sample.  Ln-material  gages, 
however,  present  their  own  class  of 
problems,  some  of  which  are  discussed 
in  Ref.  18.  Included  (for  the  carbon  gage) 
are  nonlinear  response,  response  time  of 
approximately  0.1  /us  which  increases  at 
stresses  less  than  0.5  GPa,  and  the  likeli- 
hood of  a previously  unrecorded  phenom- 
enon of  a change  in  piezoresistive  coef- 
ficient when  the  gages  are  used  in  porous 
materials  which  subject  the  gage  to  micro- 
stretching. Application  of  in-material 
gages  in  porous  samples  is  a relatively 
new  concept.  Those  interested  in  further 
details  should  consult  Refs.  3 and  18. 

The  gages  measured  wave  profiles  in 
the  form  of  voltage-time  oscillograms, 
which  were  later  converted  to  stress-time 
histories  using  a voltage-stress  calibra- 
tion procedure.  The  accuracy  of  the  wave 
profiles  measured  using  in-material  gages 
was  affected  by  the  following  sources: 

(i)  possible  tilting  of  an  impactor  relative 
to  the  target  at  the  time  of  impact, 

(ii)  nonplanarity  (hills  and  valleys)  of  the 
impact  surface  of  the  porous  specimens, 
and  (iii)  local  inhomogeneities  through 
which  the  waves  propagated. 

Most  of  the  wave  arrival  time  errors 
originating  from  the  tilt  were  corrected 
by  using  coaxial  shorting  pins  to  measure 
the  degree  of  impactor  tilt.  The  increase 
in  apparent  wave  risetime  from  the  tilt 
of  the  shock  wave  as  it  passed  over  the 


sensitive  element  of  the  gage  was  not  cor- 
rected for  in  the  data  analysis.  In  all  but 
a few  of  the  shots,  however,  the  additional 
risetime  was  less  than  20  ns  and  can  be 
ignored. 

The  second  difficulty  was  partly  re- 
moved by  applying  a thin  coat  (<2  /um)  of 
impedance-matching  epoxy  to  the  impact 
surface  to  make  a flatter  and  smoother 
surface.  Nevertheless,  the  velocities 
calculated  from  the  first  signal  arrival 
times  measured  at  different  gage  stations 
were  slightly  irregular.  To  make  a 
meaningful  comparison  between  the  exper- 
iments and  the  theoretical  calculations, 
the  experimental  wave  profiles  were 
adjusted  to  conform  to  an  average  velocity. 
For  the  sintered  materials,  the  experi- 
mental arrival  times  of  the  precursor 
were  readjusted  by  fixing  0.62  cm//us  as 
the  velocity  at  0.2  GPa  for  all  measured 
wave  profiles  except  for  the  two  thickest 
specimens  where  the  waves  were  atten- 
uated below  0.2  GPa  (in  which  case  the 
arrival  times  were  obtained  using 
1.15  cm  /us,  the  ultrasonic  velocity,  as 
the  velocity  at  the  foot  of  the  wave  fronts). 
The  value  of  0.62  cm  /us  was  chosen  by 
observing  that  the  precursor  front  main- 
tained a reasonably  constant  spreading  as 
the  precursor  traversed  the  material. 

The  "toe"  of  the  precursor  was  assumed 
to  travel  at  longitudinal  sound  speed, 

1.15  cm '/us.  Since  the  toe  was  not  easy 
to  observe  on  all  wave  profiles,  the 
stress  assumed  to  have  a constant  velocity 
was  0.2  GPa,  slightly  below  the  Hugoniot 
elastic  limit  of  the  material.  Other 
stresses  and  their  corresponding  veloc- 
ities were  tested  and  gave  approximately 
the  same  arrival  times  for  the  wave 
profiles. 
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For  the  as-sprayed  materials,  where 
no  well-defined  elastic  wave  was  detected, 
a similar  procedure  gave  less  satisfactory 
results,  Ln  this  case,  we  calculated  L' 

s 

from  the  Hugoniot  relations  and  Eq.  (lb). 

It  was  in  turn  used  as  the  compaction  wave 
velocity  at  the  midpoint  of  the  peak  stress 
level  in  the  measured  wave  profiles.  In 
general,  the  readjustment  to  the  constant- 


velocity  criterion  resulted  in  time  shifts 
of  less  than  0.09  p s , except  for  the  thickest 
(1.9- cm)  specimens,  for  which  the  time 
shifts  were  considerably  larger.  We 
noted,  however,  that  the  wave  profiles  at 
this  thickness  indicated  rarefactions  from 
the  specimen  edges  had  affected  the  wave 
forms  to  an  extent  that  made  the  time 
adjustment  of  questionable  value. 
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3.  Description  of  the  Porous  Constitutive  Model 


Our  model  is  based  on  Holt  et  al.'s 
P-o-t  model  which  is  in  turn  a rate- 
dependent  modification  of  Herrmann's 
model.  The  new  features  in  our  model 
are:  (i)  inclusion  of  deviatoric  stresses, 

(ii)  allowance  for  partial  reopening  of  the 
pores  upon  unloading,  (iii)  use  of  an  im- 
proved expression  for  the  static  compaction 
function,  and  ( iv)  use  of  a porosity- 
dependent  relaxation  time  for  pore  closure. 

3.1  FIRST  CONSTITUTIVE 
RELATION:  EOS  OF 
POROUS  MATERIALS 

Two  relations  are  used  in  the  formu- 
lation of  the  constitutive  model.  The  first 
relation  is  an  equation  which  relates  the 
pressure  (P),  volume  (V),  and  energy  (E) 
of  a porous  material  to  the  equation  of 
state  (EOS!  of  the  corresponding  fully 
compacted  solid  (hereafter  denoted  simply 
as  "solid"  for  brevity).  The  porosity 
parameter  a is  defined  by  this  relation: 

P(V,  E)  i Ps(v  <*,  Ks>,  E = Es,  (2a) 

which  can  be  rewritten  in  an  alternate  but 
equivalent  form, 

a - V V = P P,  E = E , (2b) 

s s'  ’ s 

where  the  subscript  s refers  to  the  solid 
material. 

The  first  equality  in  Eq.  (2b)  states  that 
the  porosity  parameter  a is  the  ratio  of 
the  porous  volume  to  the  solid  volume. 

The  second  equality  in  Eq.  (2b)  states  that 
in  taking  the  volume  ratio  the  porous 
volume  should  be  calculated  at  the  porous 
pressure  P which  is  smaller  than  Pg  by 
a factor  l/o.^  (See  points  C and  C"  in 


f ig.  5(a).)  This  is  a physically  reasonable 
procedure,  since  the  pressure  is  the  force 
per  unit  cross-sectional  area,  and  since 
the  cross-sectional  area  for  the  porous 
material  is  larger  by  a factor  a. 


Volume 


Volume 

(b) 


Fig.  5.  (a)  Schematic  diagram  showing 

porous  and  solid  hydrostats. 
Loading  and  unloading  paths  for 
a porous  material  are  also  indi- 
cated. (b)  Dynamic  paths  of  a 
sintered  porous  sample  are 
schematically  shown.  Note  that 
these  paths  are  not  unique,  and 
that  wave  profiles  measured  at 
points  close  to  the  impactor- 
target  interface  should  qpntain 
more  pronounced  time-dependent 
effects  than  those  obtained  from 
a transmitted  wave  experiment. 


-13- 


" 


3.2  SECOND  CONSTITUTIVE 

HEI.ATION:  EQUATION  EUR 
THE  POROSITY  PARAMETER 

The  second  relation  is  an  equation 
which  describes  a path  for  o during  a 
time-dependent  deformation  process. 

Using  the  porous  beryllium  hydrostats  in 
Eig.  2 as  a guide,  one  can  distinguish 
essentially  two  different  types  of  paths 
that  a must  follow,  namely  plastic  and 
elastic  paths.  The  terms  plastic  and 
elastic  used  here  refer  to  the  irreversible 
and  reversible  portions  of  the  pressure- 
volume  paths  under  hydrodynamic  (or 
hydrostatic)  deformation.  The  irrever- 
sible behavior  of  the  P-V  path  is  charac- 
teristic of  a porous  material  and  is 
associated  with  a microscopic  plastic 
flow  of  material  adjacent  to  the  pores, 
where  local  stress  concentration  is  very 
large.  It  must  be  distinguished  from  an 
additional  elastic-plastic  effect,  arising 
from  a macroscopic  shear  stress  and 
subsequent  yield  behavior,  which  is  taken 
into  account  as  the  deviatoric  stresses 
and  is  considered  later.  Figure  5(a)  illus- 
trates the  elastic  and  plastic  paths  under 
hydrostatic  deformation. 

For  sintered  beryllium,  loading  follows 
initially  along  an  elastic  path  (AjA')  until 
pressure  reaches  the  "yield  pressure" 

P^,,  where  the  microscopic  plastic  flow 
process  begins  to  dominate  over  the  local 
elastic  strain  contribution  to  P.  We 
assume  that  this  path  is  reversible  upon 
unloading  and  is  characterized  by  a rela- 
tively small  compressibility  on  account  of 
a large  number  of  more  spherical  (hence 
stronger)  pores  present  in  the  sintered 
specimen.  Beyond  point  A',  the  deforma- 
tion follows  a plastic  path  (A'B'C'),  which 
is  mostly  irreversible.  Along  this  por- 


tion of  the  P- \ path  the  specimen  is  more 
easily  compressed.  The  implication  of  a 
negative  curvature  at  point  A'  for  wave 
profile  measurements  is  that  such  curva- 
ture introduces  a precursor  wave  which 
precedes  the  main  shock  wave  associated 
with  the  plastic  path.  In  accordance  with 
the  static  experiments,  unloading  from 
and  reloading  back  to  the  plastic  path 
(Al  - A',  B *-  B',  and  C - C'  in 
Fig.  5(a)),  as  well  as  tensile  loading 
(P  < 0),  are  assumed  to  follow  elastic 
paths. 

For  as-sprayed  porous  beryllium,  the 
irregularly  shaped  pores  and  cracks 
present  are  more  easily  deformed  (thus 
more  easily  closed)  than  the  more 
spherically  shaped  pores  in  the  sintered 
material  (with  which  elastic  deformation 
is  associated).  Therefore,  plastic  loading 
in  the  as-sprayed  beryllium  is  assumed 
to  start  essentially  at  zero  pressure,  i.e., 
from  point  A^  in  Fig.  5(a).  Such  a P-V 
path  will  produce  a wave  profile  without 
a precursor.  Aside  from  this  difference, 
deformation  characteristics  of  the  as- 
sprayed  specimen  are  taken  to  be  quali- 
tatively similar  to  those  of  the  sintered 
specimen,  although  the  parameters  in  the 
model  are  quite  different  from  each  other. 

3.2.1  Plastic  Path 

During  dynamic  loading  processes  such 
as  those  shown  in  Fig.  5(b),  where  the 
pressure  and  volume  of  porous  materials 
lie  above  the  static  compaction  path  in 
Fig.  5(a)  (curves  A^A'B'C'  and  A'B'C'  for 
as-sprayed  and  sintered  beryllium,  re- 
spectively), and  a is  still  decreasing  with 
time,  we  assume  the  porous  materials 
follow  a plastic  path.  Along  this  path  we 
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ust-  the  rate-dependent  pore-closure  rela- 
tion  of  Holt  et  al., 

o g(P)  - t do  dt.  (3) 

The  parameter  t is  a time  constant 
describing  the  rate  of  plastic  flow  of 
material  into  pores,  which  reduces  o to 
a final  average  equilibrium  value  given 
by  the  static  compaction  function  g(P). 

To  determine  g(P),  we  eliminate  V in 
Eq.  (2)  by  using  hydrostatic  P-V  data  for 
porous  specimens  and  similar  P-V  data 
for  solid  beryllium. *‘  The  resulting  data, 
P and  g(P),  can  be  accurately  represented 
by  the  following  expression: 

g(P)  -ur  + (on  - ft  ) exp(aP  + bP“  + cP'5), 

or  U cr 

(4) 

where  ctn  is  the  initial  porosity  and  ft  is 

introduced  to  account  for  some  residual 

porosity  that  appears  to  persist  under 

7 18 

extremely  high  pressures.  ’ The  con- 
stants a,  b,  and  c are  obtained  by  fitting 
the  data  by  a least-squares  method.  Other 

mathematical  representations  used 
8 9 

earlier,  ’ namely  polynomials  in  P,  did 
not  represent  the  experimental  data  as 
well  as  Kq.  (4)  does.  It  is  perhaps  im- 
portant to  note  that  Eq.  (4)  is  derived 
from  the  hydrostatic  data;  yet,  as  we 
shall  see  later,  it  plays  a major  role  in 
reproducing  dynamic  data  such  as  wave 
profiles. 

There  is  a close  similarity  between 
Eq.  (3)  and  another  rate-dependent  relation 
derived  from  a spherical  model  calculation 
for  an  elastic-viscoplastic  material  by 
Holt  et  al.  * ' ; i.e., 

P = Pstatic(a)  - T'  d"/dt' 

t’  r 4n/[3 ft(ft  - 1)]  . (5) 

In  deriving  Eq.  (5),  the  inertial  term, 

-1 


which  is  small*  compared  to  the  viscous 
contribution  i 1 do  dt,  was  neglected. 
Equation  (5)  describes  an  approach  of  the 
dynamic  pressure  P(a,  t)  to  the  static 
value  Psta^.c<ft),  where  a is  an  independ- 
ent variable,  while  Eq.  (3)  describes  an 
approach  of  the  dynamic  ft(P,  t)  to  the 
static  value  g(P),  using  P as  an  inde- 
pendent variable.  Since  the  two  expres- 
sions represent  the  same  relaxation 
process,  the  time  constants  i and  t'  will 
likely  be  equal  or  differ  at  most  by  a 
multiplicative  constant.  Assuming  that 
this  is  true,  and  replacing  the  viscosity  rj 
in  Eq.  (5)  in  terms  of  ftQ  and  ( t at  o^), 
we  obtain  the  following  ft-dependent  t: 

T = To  (Qo  * QX)^Q{Q  - %>)•  (6) 

Again,  a ^ is  used  in  place  of  unity  to 
account  for  the  observed  behavior  of  a at 
high  pressures. 

Equation  (6)  gives  a value  for  t which 
is  equal  to  at  the  foot  (P  - 0)  of  the 
compressive  wave  front,  but  which  becomes 
larger  at  higher  stress  levels.  By 
choosing  small  enough  we  can  obtain 
calculated  wave  profiles  for  as-sprayed 
beryllium  similar  to  those  observed  in 
plate  impact  tests,  namely  a steeply  in- 
creasing portion  at  the  foot  which  may 
become  ramped  (depending  on  P)  at  the 
peak  stress  level.  It  should  be  noted  that, 
for  physically  reasonable  systems,  al- 
though r becomes  large  at  high  pressures, 
r dft/dt  in  Eq,  (3)  always  approaches  zero. 

3,2.2  Elastic  Path 

If  the  conditions  for  the  plastic  loading 
path  are  not  satisfied,  the  porous  material 
is  assumed  to  follow  an  elastic  path,  where 
the  pores  respond  instantaneously  without 


t 


any  time  delay.  Assumption  of  the  time- 
independent  elastic  path  has  been  checked 
to  be  nearly  correct  in  the  case  of  a hollow 
sphere  in  an  aluminum  matrix  “ 

To  obtain  an  equation  for  a along  an 
elastic  path,  we  first  express  the  porous 
and  solid  EOS's  in  incremental  form. 

That  is,  using  A to  denote  the  difference 
in  P (and  in  P ) at  two  successive  time 
steps  (differencing  as  used  in  a numerical 
hydrodynamic  code),  we  have 


AF  = MAp  p). 

(7a) 

AF  = K (Ap  , p ), 
s s ^s  rs 

(7b) 

where  K and  Ks  are  the  bulk  moduli  of 
the  porous  and  solid  materials,  respec- 
tively. The  density  p in  Eq.  (7a)  is 
eliminated  by  using  Eqs.  (2a)  and  (7b). 
After  lengthy  algebraic  manipulation,  we 
obtain  an  expression  for  the  rate  of  change 
of  o with  P along  the  elastic  path;  i.e., 

da/dP  = a(K  K - a)/(aP  - K ).  (8) 

s s 

In  a numerical  hydrodynamic  code,  the 

use  of  an  incremental  expression  such  as 

Eq.  (8)  rather  than  a specific  form  for  a 

is  preferred.  That  is,  during  the  initial 

compaction,  a switches  from  the  elastic 

ae  given  by  a numerical  integration  of 

Eq.  (8)  to  the  plastic  given  by  solving 

Eq.  (3)  whenever  a exceeds  a . No 
e p 

additional  parameter,  such  as  the  pres- 
sure at  the  elastic  limit  (point  A'  in 
Fig.  5(a)),  is  needed.  Also,  since  Eq.  (8) 
does  not  require  K and  K to  be  linearly 
elastic,  any  desired  forms  for  K and  Kg 
can  be  inserted  whenever  necessary. 

Next,  the  bulk  modulus  K required  in 
Eq.  (8)  and  the  shear  modulus  G required 
in  the  deviatoric  stress  must  be  specified. 


For  this  purpose,  Mackenzie's  expres- 
sions^ for  K and  G,  namely 


1 K 


Ks  4 <(4G  11  4 °l(a  - (9a) 


5G  (K  + 4G  ) 

G G - s s s o - 1 


ATK — r fJTT  - 

s s 


+ (.)( (a  - 1)1,  (9b) 


for  a linear  elastic  medium  containing 
noninteracting  spherical  pores  are 
rewritten;  i.e.. 


1 /K  = (a^/Kglt  + <1/K0><1  - t),  (10a) 

G = <1/Gg)t  + <1/G0>(1  - t), 

t (»0  - Q')/(o0  - a^),  (10b) 

where  the  subscript  0 refers  to  the  initial 
porous  state.  These  expressions,  which 
were  first  introduced  by  Seaman1 1 (with 
a v - 1),  are  obtained  directly  by  elimi- 
nating G in  Eq.  (9a)  in  terms  of  Kn  and 

o u 

and  Kg  in  Eq.  (9b)  in  terms  of  G^ 
and  o^.  The  factor  (1  - t)  in  Eq.  (10b)  is 
actually  (1  - t)a  /a  which  is  approximated 
to  (1  - t)  + 0|(a  - l)2] . 

The  use  of  Eq.  (10)  rather  than 
Mackenzie's  original  expressions  (Eq.  (9)) 
is  clearly  preferred,  since  Eq.  (10)  is 
designed  to  reproduce  the  correct  moduli 
at  two  limiting  values  of  a,  namely 
and  a The  straightforward  use  of 
Mackenzie's  expressions,  Eq.  (9),  which 
are  correct  to  0((o  - l)2)  for  K and 
O (a  - 1)  for  G,  will  yield  poor  results  for 
a material  with  a large  porosity  as  well 
as  for  materials  such  as  as-sprayed 
beryllium  whose  initial  moduli  KQ  and  GQ 
are  greatly  different  from  Kg  and  Gg, 
respectively. 


L 
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3.3  WELD  STKKS.S 

To  complete  the  formulation  the  devia- 
toru-  stress  s must  be  specified.  Then, 
the  axial  stress  (along  direction  1)  in 
a 1-1)  strain  deformation  is  given  by 

a , * P * Sj,  (11) 

where 

Sj  ^-Gln(V  YQ),  if  |sj|  < j V,  (12a) 


t ~ Y otherwise  (+  for  loading  and 

- for  unloading).  (12b) 


The  shear  modulus  G in  Eq.  (12a)  is 

calculated  from  Eq.  (10b).  Y may  be 
2 

calculated  from  Y = *■  P if  experi- 

mental data  are  available  for  P and  Oj. 
The  data  can  be  fitted  by  a functional 
form  such  as 


Y = Max  [V  Y1  + V2C 


(13) 


where  c = IhIV/Vq)  and  small  elastic 
strain  is  neglected.  We  use  ln(V/Vg)  for 
c rather  than  the  equivalent  plastic  strain 
invariant  t , which  is  usually  used  in  cal- 
culation of  Y for  solid  materials.  In 
most  elastic-plastic  theories  of  solids, 

c is  assumed  to  be  a nondecreasing 
P 

function  of  the  strain  history,  and  the 

resulting  Y becomes  also  a nondecreasing 

function  of  e , which  is  contrary  to  the 
P 

cyclic  loading  and  unloading  data  of 

5 

Schock  et  al.  mentioned  in  Sec.  2.2. 

Apart  from  the  parameter  t^,  the 
formulation  discussed  above  requires 
only  static  compression  data  and  ultra- 
sonic data  at  P = 0 to  be  measured  in  the 
porous  material,  and  all  other  parameters 
can  be  determined  from  the  properties  of 


the  solid  materials,  in  practice,  such  a 
procedure  has  not  been  successfully 
applied  beyond  1 GPa,  mainly  because 
the  1-D  strain  experiment  from  which  <; ^ 
i.s  obtained  requires  a sophisticated  tech- 
nique of  applying  both  confining  stress 
and  axial  stress  and  is  difficult  to  perform 
at  higher  stresses.  Therefore,  at  pres- 
sures greater  than  the  experimental  limit 
of  our  static  1-D  strain  apparatus  (-  1 GPa) 
the  Iiugoniot  stress  (Ojj)  rather  than 
was  used  to  calculate  Y.  In  the  case  of 
as- sprayed  beryllium  the  data  on  and 
agree  nearly  perfectly  to  1 GPa,  thus 
providing  an  experimental  support  for 
using  C7j|  in  place  of  at  higher  stresses. 
We  note  that  the  yield  strength  obtained  in 
this  manner  exceeds  the  value  (-0.2  GPa) 
corresponding  to  solid  beryllium  at  some 
partially  compacted  volumes.  Such  a 
phenomenon  would  seem  to  be  entirely 
possible,  since  the  yield  stress  can  in- 
crease through  both  the  work  hardening  of 
the  matrix  material  and  the  densification 
during  deformation.  At  some  region  of 
partial  compaction,  therefore,  the  latter 
could  give  a dominant  contribution  to  Y as 
would  seem  to  be  the  case  here. 

The  corresponding  comparison  for 
sintered  beryllium  gives  an  apparently 
inconsistent  result.  The  comparison 
shows  that  lies  above  cr  with  the 
deviation  becoming  larger  at  higher 
strains.  At  (Vq  - V)  Vq  0.04,  for 
example,  Ojj  is  about  1.2  GPa  while  is 
about  0.9  GPa.  This  disagreement  is  dif- 
ficult to  explain  at  present.  The  correc- 
tions to  Cj  or  from  the  thermal  con- 
tributions are  too  small  to  explain  the 
differences.  Also  a possible  rate- 
dependent  correction  to  stress  from 
viscous  damping  is  a non-steady-state 
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process  which  would  affect  only  the  shape 
of  the  shock  front,  not  the  amplitude. “ 1 
The  Hugoniot  stresses  were  taken  well 
behind  the  shock  front  and  would,  there- 
fore, not  be  affected  to  a significant 
extent  by  the  rate-dependent  process. 
Further  experimental  and  theoretical 
effort  would  seem  to  be  useful  to  clarify 
this  situation. 

Having  no  clear  preference  of  one 
choice  of  V over  another  for  sintered 
beryllium,  we  made  two  separate  calcu- 
lations on  wave  propagation  to  test  the 
form  of  Y.  For  the  first  calculation,  Y 
was  taken  to  be  constant  (0.1  GPa).  This 
choice  was  based  on  the  assumption  that 
overall  differences  in  the  hydrostat  and 
<7j,  which  follow  each  other  very  closely 
to  the  highest  stress  level  achieved  in  the 
static  1-D  strain  experiment  (-1  GPa), 
can  be  approximately  represented  by 
using  the  constant  yield  strength  (~0.1  GPa). 

Actually  the  static  data  in  the  region 
of  compaction  near  the  yield  point  present 
a more  complex  picture  than  the  constant 
Y assumption  used  here,  since  second 
order  effects  (e.  g. , material  anisot- 
ropy) which  are  known  to  exist  in  porous 
beryllium  but  which  were  ignored  in  our 
analysis,  become  relatively  more  signifi- 


cant here.  Numerous  attempts  were 
made  to  reconcile  the  yield  behaviors  as 
measured  in  the  different  experiments. 

In  the  end,  the  pragmatic  view  was  taken 
that  the  more  sophisticated  forms  of  the 
yield,  when  placed  in  the  hydrocode,  did 
not  produce  significant  differences  in  the 
predicted  wave  profiles  investigated  here. 
At  low  stress  levels  (<0.3  GPa),  however, 
a refined  model  may  be  needed  in  the 
future  to  account  for  the  anisotropic  effect 
explicitly. 

For  the  second  calculation,  Y was  ob- 
tained from  1.5  (a„  + P),  i.e.,  Eqs.  (11)- 
(12),  using  the  data  on  and  P.  The 
resulting  numerical  values  of  Y were 
fitted  by  a polynomial  form  of  Eq.  (13). 

A test  calculation  was  carried  out  using 
shot  geometry  identical  to  shot  BU-B  (see 
Fig.  8(d)  in  the  next  section).  The  result 
showed  that  the  constant  Y method  gen- 
erally gave  much  better  agreement  with 
experiment  than  Y obtained  from  a fit  to 
the  Ojj-Pdata.  In  the  latter  case,  the 
compressive  wave  front  had  a faster  rise- 
time and  was  traveling  so  fast  that  the 
experimentally  observed  attenuation  did 
not  take  place  in  the  computer  calculation. 
The  calculations  shown  in  Sec.  5 were 
made  using  Y = 0.1  GPa. 
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4.  Computational  Procedure 


Solutions  of  Eqs.  (2a)  and  (3)  give 
o ( »)  applicable  to  the  plastic  path, 

while  Eqs.  (2a)  and  (8)  give  the  elastic 
solutions  Since  both  pairs  of  equations 
are  highly  complex,  we  solve  them  iter- 
atively in  the  hydrodynamic  code  KO.^ 
Iterative  schemes  for  both  cases  are 
similar,  however,  so  a single  scheme 
was  adopted  to  take  care  of  both  cases. ^ 

We  introduce  a function  F (P), 

a 


a , if  a < 
e’  e 


N . 

or  a < a 


P 

or  P t 0,  (14a) 


F ( P)  - x 1 , if  a or  a 


a p,  otherwise. 


1, 


(14b) 


(14c) 


where  the  superscript  N represents  the 

present  time  step  in  a finite-difference 

scheme  used  in  the  code,  and  a , a , 

e g 

and  Op  are  defined  by 


Q(i  - aN  + (P  - PN)(do/dP), 


N 


(do/dP)  = Eq.  (8)  at  a = a , (15a) 
g(PN+2),  PN+5  - (P  + PN)/2,  (15b) 

o-p  = (c/g  - o^/2  + o^T/At)/(0.5  + T/At). 


(15c) 

That  is,  Fff(P)  is  set  to  if  (i)  a lies 
below  the  static  compaction  path  a or 
(ii)  the  dynamic  value  op  at  the  next  time 
step  starts  to  decrease,  or  (iii)  the  zone 
under  consideration  is  in  tension. 

We  also  define  the  second  function 
Fp(»),  which  is  Eq.  (2a)  evaluated  at  the 
(N  +•  l)th  time  step,  except  that  the  correct 
a is  yet  to  be  determined;  i.e., 

Fp(o)  = i Ps(VN  + 1/».  EN+1).  (16) 


The  basic  idea  behind  the  scheme  is  to 
enclose  a desired  solution  by  a two- 
dimensional  "box"  with  its  sidelengths 
a"  - o'  (>0)  and  P"  - P'  (>0),  and  to 
reduce  the  size  of  the  box  (holding  the 
solution)  in  each  iterative  cycle  until  the 
solution  has  converged  within  a desired 
accuracy.  Specifically,  this  is  done  as 
follows; 

(1)  Initial  box  size: 

»’  MinjV*,  Fa(Fp(nN))], 

a"  Max[oN,  Fa(Fp(nN))J, 

P'  = Min[pN,  Fp(Fo(PN))], 

P"  = Max[pN,  Fp(F^(PN))]  . 

(2)  Iterative  loop: 

(i)  = (a1  + a")/2, 

PA  = (P>  + P")/2. 

(ii)  o-iteration  to  obtain  new  o'  and  a": 
«B  = Min[V  FjFp(«A))], 

«c  = Max[aA,  Fa(nA)]  , 

a'  = Max  Ja',  cigj, 
a " = Min  jo", 

(iii)  P- iteration  to  obtain  new  P'and  P' 

PB.M,n[pA,  Fp(F„<PA>)] . 

pc  ■ M“  [''a-  Fp(F„«pa')]. 

P'  = Max[p',  Pg], 

Pn  = Min  [P",  Pc], 

(3)  Convergence  criterion: 

The  convergence  test  for  | P'  | 

< 10'7  GPa  is 


19- 


k 


f- 


H53 


i 


i 


} 


jC,f  - ~ ,C/- ■ j f ! I'1  0.005; 

the  convergence  test  for  | P' | > 10  1 GPa 
is 


o'  - o" 

1 + 

1 p"  - P'  1 

o' 

1 

1 p.  1 

[f  the  inequality  holds,  we  take  u (o' 
r fir")  2 and  P (P1  ♦ P")  2 as  our  final 
solutions;  otherwise,  the  program  returns 
to  the  next  iterative  loop. 

Note  that  the  above  scheme  performs 
essentially  two  separate  iterations,  one 
on  o and  another  on  P,  rather  than  the 
usual  way  which  relies  on  a single  itera- 
tion loop  on  a,  since  P is  the  dependent 
variable  of  a.  This  precaution  was  nec- 
essary, for  P can  change  enormously  for 
a slight  change  in  a,  especially  near  the 


region  of  complete  compaction.  (See 
Kef.  14  for  further  details.) 

Summarized  in  Table  2 are  the  exact 
numerical  values  that  were  used  in  the 
wave-profile  calculations  for  porous 
beryllium.  A constant  r was  used  for 
sintered  beryllium,  since  it  gave  a slightly 
better  wave  profile  than  the  one  obtained 
from  Eq.  (6).  The  difference  is  minor, 
however.  The  values  of  and  Gq  are 
those  obtained  from  ultrasonic  measure- 
ments. The  expressions  for  g(P)  and  Y 
are  applicable  to  pressures  less  than 
4 GPa.  Application  of  these  formulae  to 
higher  stresses  is  subject  to  uncertain 
error. 

Calculation  of  wave  propagation  also 
requires  knowledge  of  accurate  EOS  data 


Table  2.  Numerical  values  used  in  the  wave  profile  calculations  for  porous  beryllium. 

Static  compaction  function  (Eq.  (4)) 

o = 1.02,  P < 4 GPa. 
oc 

Sintered  He;  a =0.1212  (GPa)'1,  b = -0.7994  (GPa)'2,  c 0.1361  (GPa)'1. 

As-sprayed  Be:  a = -0.4822  (GPa)'1,  b = -0.4355  (GPa)'2,  c = 0.6975  (GPa)'2. 

Relaxation  time  for  pore  closure 
Sintered  Be:  r = 0.04  ps  (constant  t). 

As-sprayed  Be:  Tq  = 0.006  ps  (porosity-dependent  t from  Eq.  (6)). 

Bulk  and  shear  moduli  (Eqs.  (10a),  (10b)) 

K = 122  GPa,  C.  = 155.6  GPa. 

Sintered  Be:  KQ  = 74  GPa,  GQ  = 108.9  GPa. 

As-sprayed  Be:  Kq  - 19.4  GPa,  Gq  0. 

Yield  stress  (Eq.  (13)) 

Sintered  Be:  Y = 0.1  GPa. 

As-sprayed  Be:  Y = Max  (0,  -0.21  t 15.75e  - 54c2),  c < 0.146, 

= 0.146,  e 0.146. 


r 


for  the  impactor,  gage,  buffer,  and  back- 
ing materials  as  well  as  the  EOS  of  solid 
beryllium.  Analytical  expressions  used 
to  represent  these  materials  are  given  in 
Table  3.  The  expressions  for  fused  quart/ 

and  PMMA  are  obtained  by  fitting  Barker 
93 

and  Ilollenbach's  data-  with  an  additional 

24  95 

small  energy-dependent  correction.  * 

The  expression  for  x-cut  quartz  was  ob- 

2 (i 

tained  by  fitting  Graham's  data.  The 
Grilneisen  7 of  x-cut  quartz  is  small  and 
is  not  well  known.  Therefore  no  energy- 
dependent  correction  was  used  for  this 
material.  The  EOS  of  the  carbon-Kapton 


gage  was  generated  by  a least-squares  fit 
routine  to  data  on  a (GPa)  vs  U„  (cm/jis) 

(o  0.22385  Up  + 0.0470  Up2).  The  data 
was  obtained  in  separate  experiments  con- 
ducted in  a fused  quartz  target  with  carbon- 

1 8 

Kapton  in-material  gages. 

Finally,  we  use  the  Grilneisen  form  of 
EOS  for  solid  beryllium  given  in  Ref.  27. 
This  expression  agrees  closely  with  a 
similar  expression  in  Ref.  17  for  a small 
E.  Since  the  energy-dependent  correction 
is  very  small  for  plate  impact  tests,  no  appar- 
ent difference  occurs  in  our  calculation  by 
using  either  this  or  the  other  expression. 


Table  3. 


Expressions  for  the  equations  of  state  (EOS)  of  the  carbon-Kapton  gage,  fused 
quartz,  x-cut  quartz,  PMMA,  and  beryllium.  P = pressure  (GPa)  and  E 
= internal  energy  per  original  volume  (GPa). 


Material 


EOS  (p  p psQ-  1) 


ps0  (g  cm3) 


Carbon-Kapton 
Fused  quartz 
X-cut  quartz 
PMMA 


P = 4.310  p + 1.904  p2  + 294.3  p3 
P = 77.56  p - 480.7  p2  + 3074  p 3 + 7.52  E 
P = 86.66  p + 45.87  p2 
P = 10.09  p + 6.272  p2  + 85  E 

p = 118,4  p (1  + 0,42  p - 0.08  p2) 

(1  - 0.124  p)2 


1.33 

2.201 

2.65 

1.184 


Be  ryllium 


+ (116  + 16  p)  E 


1.85 


5.  Comparison  Between  Predicted  and 
Experimental  Wave  Profiles 


Experimental  wave  profiles  which  will 
be  used  for  comparison  were  obtained  at 
two  stress  levels,  approximately  O.G  and 
1.7  GPa.  Profiles  were  obtained  by  im- 
pacting a PM  .VIA  plate  onto  either  sintered 
or  as-sprayed  porous  beryllium  having 
carbon- Kapton  gages  embedded  at  as  many 
as  six  different  depths  within  the  material. 
Thus  a single  experiment  took  data  on  wave 
profiles  at  several  different  levels.  The 
gages  were  located  as  deep  as  1.9  cm  from 
the  impact  surface  and  some  measurement 
times  extended  up  to  4 ps.  The  resulting 
21  wave  profiles  — both  unattenuated  and 
attenuated  — represent  large  variations  in 
the  physical  parameters  and  enabled  us  to 
extensively  test  the  model.  Experimental 
parameters  (shot  names,  thicknesses  of 
impactors,  target  specimens,  and  gages, 
measured  peak  stresses,  durations  of  the 
measurements,  etc.)  are  summarized  in 


Tables  4 and  5 for  the  sintered  and  as- 
sprayed  porous  specimens.  In  the  compu- 
tation the  thicknesses  of  the  backing  mate- 
rials (which  were  placed  behind  the  gages 
to  carry  off  the  wave  without  reflection) 
are  reduced  from  the  original  (experi- 
mental) thicknesses  in  some  cases 
(Tables  4 and  5)  since  this  saved  the 
computing  time  without  affecting  the  cal- 
culated stress  histories. 

Measured  and  computed  wave  profiles 
are  compared  in  Figs.  6,  7,  and  8 for  the 
sintered  material  and  in  Figs.  9,  10,  and 
11  for  the  as-sprayed  material.  Note 
that  abridged  notations  such  as  shot  BT-A 
and  shot  BT-B  are  used  here  to  indicate 
that  the  stress  profiles  were  recorded  at 
levels  A and  B in  shot  BT.  Agreement 
between  the  experimental  wave  profiles 
and  those  predicted  from  the  present 
model  is  good  both  quantitatively  and 


(a)  Time  (nsec)  (b)  Time  (nsec) 

Fig.  6.  Computed  and  measured  full  wave  profiles  for  the  sintered  porous  Be  specimen 
at  (a)  high  and  (b)  low  stress  levels.  See  Table  4 for  experimental  parameters 
used  in  these  shots. 
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Fig.  7.  Attenuated  wave  profiles  of  the  high-stress  shot  (shot  BS)  for  sintered  porous 
Be:  (a)  experimental  gage  records,  (b)  impactor-target  geometry  (see  also 
Table  4),  (c)-(d)  comparison  of  measured  and  computed  wave  profiles  at  dif- 
ferent target  thicknesses.  (See  text  for  possible  corrections  to  shot  BS-B.) 


| 

qualitatively.  We  give  below  more  de- 
tailed comparisons  of  the  compressive 
and  release  portions  of  these  wave  pro- 
files and  discuss  several  interesting 
features,  including  the  effect  of  gage 
interaction  on  the  wave  profiles. 


5.1  COMPRESSIONAL  W AVES 

For  sintered  porous  beryllium  (Figs.  7 
and  8),  an  elastic  precursor  about  0.4  GPa 
high,  which  corresponds  to  the  "shoulder" 
in  the  hydrostat  in  Fig.  2,  precedes  the 
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Fig.  8.  T ull  and  attenuated  wave  profiles  of  the  low-stress  shot  (shot  BU)  for  sintered 
porous  Be:  (a)  experimental  gage  records,  (b)  impactor-target  geometry  (see 
also  Table  4),  (c)-(f)  comparison  of  measured  and  computed  wave  profiles  at 
different  target  thicknesses.  (See  text  for  possible  corrections  to  shot  BU-D.) 


main  plastic  wave.  At  the  foot  of  the  pre- 
cursor the  velocity  is  close  to  the  longi- 
tudinal sound  speed,  but  it  becomes  slower 
at  higher  stresses.  The  value  0.62  cm/ps 
at  P = 0.2  GPa  chosen  for  the  calculation 
of  the  first  arrival  times  agrees  reason- 


ably well  with  the  velocity  of  the  precursor 
at  this  level. 

Figures  12  and  13  give  comparisons  of 
the  development  of  the  non- steady -state 
shock  wave  and  the  precursor.  These 
transmitted  compressive  profiles  were 


-26- 


0 2 1 1 1 1 1 1 1 1 1 1 I I 1 I l l l l l l 

'0  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6  1.8  0 0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6  1.8 


(a)  Time  (usee)  (b)  Time  (usee) 

Fig.  9.  Computed  and  measured  full  wave  profiles  for  as-sprayed  porous  Be  at  (a)  high 
and  (b)  low  stress  levels.  See  Table  5 for  experimental  parameters  used  in 
these  shots. 

measured  at  four  different  thicknesses  by  A faster-rising  shock  front  and  the 

x-cut  quartz  gages.  The  impactor-target  lack  of  an  elastic  precursor  are  two  chief 
arrangement  of  Fig.  3(b)  was  used  here.  distinguishing  features  of  compressive 

We  note  that  agreement  between  the  pre-  profiles  in  as-sprayed  porous  beryllium, 

dieted  and  experimental  wave  profiles  is  Good  agreement  along  the  compressive 

satisfactory.  portion  of  the  wave  has  largely  resulted 
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Fig.  10.  Attenuated  wave  profiles  of  the  high-stress  shot  (shot  BT)  for  as-sprayed 
porous  Be:  (a)  experimental  gage  records,  (b)  impactor-target  geometry 
(see  also  Table  5),  (c)-(h)  comparison  of  measured  and  computed  wave  pro- 
files at  different  target  thicknesses.  (See  text  for  possible  corrections  to 
shots  BT-A  and  BT-B.) 


from  the  use  of  a new  static  compaction 
function  having  an  exponential  form  rather 
than  earlier  expressions^'  **  which  use 
polynomials  of  (Pc  - P)  '(P  - P ) or 

(pc  - p)/(p  - p ).  (The  subscript  c indi- 


cates full  compaction,  and  e indicates 
elastic  limit,  i.e.,  point  A'  in  Fig.  5(a).) 
For  materials  which  stay  porous  over  a 
relatively  large  range  of  pressure,  a 
single  polynomial  expression  (even  as 


■R 


IL 


(g)  Time  (usee) 

Fig.  10. 

high  as  fifth-order)  is  not  adequate  to 
represent  the  static  compaction  curve  over 
the  pressure  range  of  interest. 

We  also  note  from  Figs.  6-11  that  the 
arrival  times  of  the  shocks  (and  their 
precursors  in  the  case  of  sintered  speci- 
mens) agree  satisfactorily.  The  relaxa- 
tion time  r = 0.04  ps  matches  the  observed 
risetimes  of  the  shocks  in  the  case  of  the 
sintered  specimens,  while  the  porosity- 


(Continued) 

dependent  r from  Eq.  (6)  adequately 
describes  the  risetimes  of  the  shocks  in 
the  as-sprayed  specimens.  In  the  latter 
case,  it  is  noteworthy  that  t changes 
from  0.060  ps  at  the  foot  of  the  shocks  to 
0.015  ps  at  1.7  GPa.  If,  instead  of  Eq.  (6), 
a constant  value  (0.015  ps)  had  been  used 
for  t,  the  foot  of  the  shocks  would  be 
traveling  too  fast,  resulting  in  much  longer 
risetimes  than  those  in  the  observed  profiles. 
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Fig.  11.  Attenuated  wave  profiles  of  the  low-stress  shot  (shot  BR)  for  as-sprayed 
porous  Be:  (a)  experimental  gage  records,  (b)  impactor-target  geometry 
(see  also  Table  5),  (c)-(g)  comparison  of  measured  and  computed  wave  pro- 
files. (Sec  text  for  possible  corrections  to  shots  BR-A  and  BR-C.) 


The  largest  deviation  occurs  in  the 
calculation  of  the  peak  stresses,  with  the 
experimental  data  lying  about  10%  or  less 
below  computer  predictions  (except  for 
shot  BT-A  in  Fig.  10  and  shots  BR-A  and 
BR-C  in  Fig.  11).  Fither  experimental 
uncertainties  or  approximations  in  the 
model  or  a combination  of  the  two  could 
account  for  these  relatively  small  dif- 
ferences. Experimental  profiles  of  shots 


BT-A  and  BR-A  are  anomalous  in  that  the 
peak  stresses  recorded  in  the  gages  lie 
below  the  peak  stresses  recorded  deeper 
within  the  material.  (Compare  the  peak 
stress  levels  between  shots  BT-A  and 
BT-B  in  Fig.  10  or  shots  BR-T  and  BR-B 
in  Fig.  11.)  Microstretching  of  the  gages 
or  macrostretching  of  the  gage  leads  may 
account  for  an  apparent  lower  stress  at 
this  surface. 
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Fig.  11.  (Continued) 


5.2  GAGE  EFFECT 

It  is  difficult  to  carry  out  a meaningful 
analysis  of  wave  profiles  without  consider- 
ing the  presence  of  the  in-material  (carbon- 
Kapton)  gage.  Although  the  gage  itself  is 
very  thin  (<0.01  cm),  the  discontinuity  at 


the  plane  of  the  gage  alters  the  wave  pro- 
files in  several  important  ways. 

First,  when  a relatively  sharp  com- 
paction front  encounters  the  gage  (such  as 
observed  in  shots  BT-A,  BR-A,  and  BIl-C), 
the  gage  response  is  not  instantaneous, 
but  because  the  wave  has  to  travel  through 


■ 


X (cm)  f . 

\ L T 

A.  0.1008  J — Theory 

B.  0.7137  // 


are  likely  causes  of  the  slow  risetimes  in 
the  experimental  wave  profiles  of  shots 
BT-  B (Fig.  10(d))  and  BR-C  (Fig.  11(e)). 

It  is  important  to  note  that  a ramped 
character  in  the  compressive  portion  of 
the  computed  full  wave  profiles  (Figs.  6, 
8(c),  9,  10(d))  was  formed  not  just  by  the 
dynamic  (t-dependent)  effect  of  the  pore 
closure  alone,  but,  as  described  above, 
nonzero  gage  ring-up  time  would  have  also 
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Fig.  12.  Comparison  of  measured  and 

computed  transmitted  compres- 
sive wave  profiles  for  sintered 
porous  Be.  Impactor-target 
geometry  used  in  this  shot 
(shot  BG)  is  also  shown.  (See 
Table  4 for  details.) 


a finite  thickness  of  the  gage  it  takes  a 

nonzero  "ring-up"  time  to  bring  its  stress 

to  that  of  the  surrounding  material.  The 

computer  calculations  give  the  ring-up 

time  as  about  0.15  qs  or  less,  which  is 

consistent  with  most  of  the  experimental 

profiles.  However,  the  experimental 

ring-up  time  may  become  even  larger  than 

the  computer  prediction  under  some  cir- 
1 8 

cumstances.  For  example,  the  carbon- 
Kapton  gages  were  found  to  respond  par- 
ticularly slowly  (0.5  ps)  at  stress  0.5  GPa 
in  a test  series  in  (solid)  fused  quartz. 
Ring-up  time  and  such  delayed  response 
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(shot  BH)  is  also  shown.  (See 
Table  4 for  details.) 
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contributed  to  shape  such  profiles, 
remark  is  based  on  our  experience  that 
the  change  of  r alone,  no  matter  how 
large,  is  not  sufficient  to  reproduce  ex- 
perimental compressive  profiles.  For 
example,  a choice  of  t which  gave  a good 
agreement  at  the  upper  portion  of  the  wave 
fronts  gave  a poor  wave  shape  at  the  foot 
of  the  wave  fronts. 

It  is  interesting  to  note  that  the  com- 
puted wave  profiles  can  reproduce  quali- 
tatively small  bumps  and  humps  atop  the 
peak  stress  levels  in  F igs.  G and  9.  Care- 
ful examination  of  the  wave  propagation  at 
discrete  time  intervals  has  shown  that 
these  minor  stress  changes  resulted  from 
an  impedance  mismatch  at  the  gage- 
specimen  interface.  A small,  early- 
arriving  portion  of  the  shock  is  reflected 
back  as  a release  wave  as  the  gage  is 
ringing  up  to  a peak  stress  level.  The 
newdy  created  release  wave  interacts  with 
the  oncoming  main  shock,  reducing  the 
latter's  stress.  For  a thin  porous  Be 
specimen,  the  release  wave  can  even 
reach  the  impact  surface  and,  upon  reflec- 
tion, create  small  compressive  waves  as 
seen  in  Figs.  6 and  9.  Disagreement 
between  the  measured  and  the  predicted 
profiles  for  waves  propagated  over  long 
distances,  e.g.  shots  BS-B  (Fig.  7(d))  and 
BU-D  (Fig.  8(f)),  originates  partly  from 
the  use  of  0.1  GPa  to  describe  an  "average" 
behavior  of  the  yield  stress  over  the  entire 
stress  range  considered  here.  This  over- 
simplification most  affected  those  cases 
where  stress  pulses  propagated  over  con- 
siderable distances  (~1.9  cm),  and  under- 
went an  appreciable  amount  of  attenuation 
from  the  peak  levels  (1.4  and  0.7  GPa, 
respectively).  The  observed  disagreement 
(-0.1  GPa)  could  have  been  reduced  with 
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more  refinement  of  the  model  parameters 
at  low  stress  levels  below  about  0.2  GPa. 

A more  important  source  of  the  disagree- 
ment, however,  occurs  from  release  waves 
that  were  created  at  the  circumference  of 
the  target  assembly.  These  waves  had 
enough  time  to  re  ich  the  gage  locations 
and  interfered  destructively  with  the  on- 
coming compressive  waves.  Analysis  of 
shots  BS-B  and  Bl'-D  with  a one- 
dimensional  hydrodynamic  code  is  there- 
fore not  well  justified. 

In  Fig.  14,  two  computed  wave  pro- 
files corresponding  to  shot  BR-D 
(Fig.  11(f))  are  superimposed  along 
with  the  experimental  wave  profile.  One 
of  the  computed  profiles  was  obtained 
including  the  gage  in  the  calculation, 
while  the  other  was  taken  at  the  same 
position  but  without  the  gage.  The 
importance  of  including  the  gage  in  a 
numerical  wave- propagation  calculation 
shows  up  rather  clearly  in  this  ex- 
ample. First,  as  might  have  been 
expected  from  the  finite  ring-up  time 
discussed  earlier,  the  presence  of  the 
gage  reduces  the  peak  stress  of  an 
attenuated  wave,  by  about  0.03  GPa  in 
this  case.  The  result  will  be  substan- 
tially increased  predicted  attenuation. 
Second,  the  gage  tends  to  smooth  out 
the  rather  pronounced  release  charac- 
teristics that  are  calculated  for  the  case 
without  the  gage.  This  release  behavior 
gives  a false  impression  that  the  mater- 
ial seems  to  be  "suddenly"  slowing  down 
the  release  wave.  As  is  clear  from 
Fig.  14,  this  is,  of  course,  not  true; 
rather,  the  behavior  is  due  to  a time 
delay  for  the  gage  to  "ring  down," 
a process  opposite  to  the  ring-up 
process  discussed  earlier. 


Stress  (GPa) 


06 


Time  [/jl sec) 

I ig.  14.  Comparison  of  measured  (dashed  line)  and  two  computed  wave  profiles  for 
shot  BR-D  (sintered  porous  Be):  solid  line  is  computed  profile  with  the 
carbon- Kapton  gage,  and  dotted  line  is  computed  profile  without  the  gage. 
(See  Table  5 for  experimental  parameters.) 


5.3  RELEASE  WAVE 

The  attractiveness  of  the  model  shows 
up  more  distinctly  along  the  release  wave 
path.  In  contrast  to  the  earlier  (hydro- 
dynamic)  models,  the  present  model  is  an 


elastic-plastic  model.  Hence,  the  experi- 
mental release  wave  velocity,  which  lies 
close  to  the  longitudinal  sound  speed,  can 
be  successfully  reproduced  in  the  model 
calculation  (Figs.  6-11).  Both  the  agree- 
ment in  the  arrival  times  of  the  release 


waves  and  the  general  shapes  of  attenuated 
wave  profiles  are  satisfactory. 

There  is  a tendency  for  the  computed 
release  wave  profiles,  after  arriving  at 
correct  times  at  the  gages,  to  release 
slightly  too  fast  in  later  portions  of  the 
release  profile.  This  can  be  seen  in 
Figs.  6,  8(e),  9,  and  10(d).  This  dis- 
crepancy may  originate  from  (i)  inaccurate 
release  behavior  in  the  model,  (ii)  erro- 
neous EOS  parameters  for  the  gage,  and 
(iii)  hysteresis  in  the  gage.  Inaccuracy 
in  the  model  probably  did  not  introduce  a 
significant  error,  since  arrival  times  of 
the  release  wave  were  accurately  pre- 
dicted from  the  model.  On  the  other  hand, 
we  observed  earlier  that  the  release  pro- 
files of  shot  BR-D  (Fig.  14)  are  affected 
sensitively  by  the  EOS  of  the  carbon- 
Kapton  gage.  This  is  also  very  likely 
true  in  the  cases  considered  here.  There- 
fore, a refined  EOS  of  the  gage  would 
likely  bring  the  measured  and  the  calcu- 
lated release  profiles  closer.  Another 
source  for  the  observed  deviation  is  the 
possibility  of  gage  hysteresis,  resulting 
in  different  resistance-stress  calibration 

scales  along  the  compressive  and  the 

1 8 

release  portions  of  the  stress  path. 

Figures  10  and  11  exhibit  a complex 
release  behavior  for  the  as-sprayed  spec- 
imen. Examination  of  the  corresponding 
model  calculations  shows  that  the  release 
of  a propagating  wave  occurs,  crudely 
speaking,  in  a stairstep  manner;  i.e., 
elastic  release  of  an  initial  peak  stress 
to  a level  calculated  from  the  von  Mises 
yield  condition,  further  elastic  release  of 
the  second  stress  level  to  the  third  stress 
level  specified  by  the  yield  condition,  and 
so  on.  Actually,  the  computed  wave  pro- 
files look  more  complicated  than  this. 


since  the  newly  created  waves  are  travel- 
ing at  different  speeds,  and  wave  inter- 
actions occur  among  such  waves  as  well 
as  among  these  waves  and  waves  reflected 
from  the  impact  surface  and  from  the  dis- 
continuities at  the  gages.  It  is  interesting 
to  note  that  this  behavior  observed  for  the 
as-sprayed  specimen  does  not  apparently 
occur  (or  is  very  small)  in  the  observed 
and  computed  profiles  of  the  sintered 
specimen  (Figs.  7(c),  8(d),  and  8(e)). 

Since  the  yield  stress  (0.1  GPa)  used  for 
the  sintered  specimen  is  relatively  small 
over  most  stress  ranges  considered  here, 
and  the  release  velocities  at  different 
stress  levels  are  nearly  equal  (~  1.2  cra'ys), 
both  the  widths  and  the  amplitudes  of 
elastic  waves  which  might  be  generated 
along  the  release  portion  of  the  wave  are 
small  and  are  smoothed  out  very  rapidly 
by  the  wave  interactions. 

5.4  OTHER  REMARKS 

Rather  than  examining  a stress  history 
recorded  at  a fixed  carbon-Kapton  gage 
station,  various  processes  discussed 
above  can  be  seen  better  by  viewing  a 
series  of  "snapshots"  of  the  wave  profiles 
taken  at  different  instants  of  time. 

Figure  15  gives  these  snapshots  corre- 
sponding to  shot  BT-E  (Fig.  10(g)).  Both 
the  elastic-plastic  character  of  the  release 
wave  and  the  effect  of  the  carbon-Kapton 
gage  can  be  seen  clearly  here.  Separation 
of  specimen  and  impactor  occurs  at  the 
impact  surface  when  the  release  wave 
from  the  back  surface  of  the  PMMA  im- 
pactor starts  to  propagate  inside  the  porous 
specimen.  The  computer-generated  wave 
profile  is  compared  with  the  corresponding 
experimental  profile  in  Fig.  15(f). 
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Fig.  15.  Snapshot  of  stress  waves  generated  from  impacting  a PMMA  plate  onto  an 
as-sprayed  porous  beryllium  target  at  impact  velocity  of  0.0750  cm  /us. 
Stress  history  obtained  at  carbon- Kapton  gage  is  shown  in  (f).  Dashed  line 
is  experimental  profile.  This  shot  corresponds  to  shot  BT-E.  (See  Table  5 
and  Fig.  10(g).) 
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Similar  snapshots  for  the  porosity 
parameter  a are  shown  in  Fig.  16.  We 
note  a large  plastic  change  in  a after  the 
passage  of  the  wave  tram.  As  shown  in 
Fig.  10(f),  this  change  occurs  because 
the  loading  and  unloading  paths  are  com- 
pletely different.  Nevertheless,  an 
appreciable  amount  of  pore  space  is  re- 
opened. This  is  shown  in  F ig.  16(c). 

It  is  easy  to  show  this  from  Eqs.  (8) 
and  (10a).  Integrating  Eq.  (8)  under  the 
assumption  of  K_.  » a P (which  is  true  in 
this  case  since  Kg  = 155.6  GPa  and 
P < 3 GPa)  gives 

O U exPl A(P,  - P))j, 

(17) 

where  P and  are  the  values  of  P and 
a at  the  initial  release  point,  and 

A <1  K0  -aQ  Ks)o.x/(a0  -aj.  (18) 

For  example,  release  from  P = 1.2  GPa 
under  static  loading  to  zero  pressure 
changes  a by  an  insignificant  amount 
(from  1.069  to  1.071)  for  the  sintered 
specimen  (o^  = 1.1254)  but  by  about  2% 
(from  1.068  to  1.090)  in  the  case  of  the 
as-sprayed  specimen  (qq  = 1.163). 

The  calculated  impact  surface  record 
on  shot  BT-A  (F'ig.  10(c))  shows  an  over- 
shoot at  the  shock  front.  Although  a 
similar  structure  was  not  found  in  the 
corresponding  experimental  record  (pos- 
sibly because  of  the  gage  ring-up  effect 
discussed  earlier),  there  are  other  ex- 
perimental records  that  show  signs  of 
such  an  overshoot.  Two  such  experimental 
records  are  shown  in  F'ig.  17,  where  a 
sintered  Be  impactor  collides  with  a buf- 
fered x-cut  quartz  target,  with  stress 
records  taken  at  0.32  cm  from  the  impact 


surface.  Unfortunately  the  tests  shown 
here  had  both  tilt  and  electrical  breakdown. 
Figure  17(a)  corresponds  to  our  model 
calculation  of  shot  BE  (Fig.  17(c))  using 
several  different  values  of  t.  These  are 
profiles  computed  not  at  the  gage  location 
but  at  the  impact  surface.  Other  addi- 
tional model  calculations  have  shown 
that  this  is  a surface  phenomenon;  i.e., 
the  overshoot  can  be  observed  only  at 
or  near  the  impact  surface  and,  at 
even  small  distances  into  the  material 
(<0.1  cm),  attenuation  has  removed  the 
overshoot  from  the  wave  front.  It  is 
easy  to  see  from  the  dynamic  path  shown 
in  Fig.  5(b)  that  such  a surface  phenom- 
enon can  be  expected  from  the  present 
model.  The  stress  at  the  impact  surface 
overshoots  because  the  stress  pulse  is 
sharply  rising  and  the  risetime  of  the 
shock  is  much  smaller  than  the  response 
time  (t)  of  the  pore  closure.  It  is  also 
noteworthy  that  the  overshoot  is  higher  at 
higher  stress  levels  (shot  BT-A  in 
Fig.  10(c)  and  shot  BE  in  Fig.  17)  and 
may  even  disappear  at  sufficiently  low 
stresses  (shot  BR-A  in  Fig.  11(c)).  It 
should  be  emphasized  here  that  the  ex- 
perimental data  presented  in  Fig.  17  are 
relatively  poor  in  quality,  so  that  any 
conclusion  based  on  them  should  be  con- 
sidered to  be  tentative  until  more  accurate 
data  become  available.  A precise  meas- 
urement of  wave  profiles  at  the  impact 
surface  would  be  useful,  since  the  decay 
time  of  the  overshoot  can  be  related  to 
t (F’ig.  17(a)).  For  such  experiments, 
the  porous  Be  specimens  must  be  care- 
fully prepared  so  that  the  surface  layers 
of  the  specimens  are  exceptionally  flat 
and  smooth,  since  the  overshoot  is  of 
such  short  duration  (<0.04  ps)  and  tilt  or 
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Fig.  16.  Snapshots  of  the  porosity  parameter  for  the  same  impactor-target  geometry 
as  used  in  Fig.  15;  (f)  shows  stress-  and  pressure-vs-strain  paths  taken 
during  the  passage  of  the  wave  train  by  Lagrangian  element  which  was  initially 
at  0.139  cm  from  the  impact  surface. 
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Fig.  17.  Stress  histories  of  sintered  porous  Be:  (a)  the  computed  profile  (shot  BE)  at 
the  impact  surface  with  different  values  of  t,  (b)  impactor-target  geometry 
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at  the  same  position  as  (c)  but  at  a slower  impactor  speed.  A = overshoot, 

B = time  of  impact,  C = stress-induced  electrical  breakdown,  D reflection 
of  wave  from  back  of  gage,  and  E = impactor  tilt  effect. 


electrical  noise  can  easily  mask  the  effect. 
The  plate-impact  geometry  of  Fig.  17(b) 
would  be  adequate  for  this  purpose,  but  a 
much  thinner  buffer  than  the  ones  used  in 


the  present  tests  would  be  desirable  so 
that  any  possible  attenuation  of  the  over- 
shoot during  the  propagation  through  the 
buffer  can  be  minimized. 
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6.  Conclusions 


In  summary,  we  have  shown  that  the 
present  model  can  describe  wave  propa- 
gation in  both  sintered  and  as-sprayed 
porous  beryllium  layers  within  the  neces- 
sary accuracy.  Both  the  compressive 
characteristics  (first-wave  arrival  time, 
risetime  of  shock,  and  peak  stress)  and 
the  release  characteristics  (arrival  time 
of  release  wave  and  general  attenuated 
wave  shape)  of  shocks  have  been  repro- 
duced satisfactorily  for  different  stresses 
and  pulse  durations  as  well  as  for  dif- 
ferent thicknesses  of  the  porous  specimens. 
Slight  deviations  between  the  experimental 
and  computed  wave  profiles  lie  mostly 
within  the  combined  errors  of  the  model 
and  the  experimental  data.  We  have  also 
shown  that  the  elastic-plastic  effect  plays 
a major  role  in  shaping  the  wave  profiles, 
and  that  release  wave  profiles  recorded 
at  the  carbon-Kapton  gage  depend  signifi- 
cantly on  the  EOS  of  the  gage.  The  model 
also  predicts  an  interesting  effect  that 
impact  surface  measurements  could  pro- 
duce, i.e.,  an  overshoot  in  the  compres- 
sive wave  profile  beyond  a certain  stress 
level. 

Important  assumptions  in  the  present 
model  are:  (i)  isotropy  of  the  medium, 

(ii)  porosity  change  directly  attributable 
to  the  pressure,  (iii)  separability  of  the 


pressure  and  the  deviatoric  stress  as 
expressed  by  Eq.  (11),  and  (iv)  a negli- 
gible Bauschinger  effect.  While  the 
isotropy  assumption  (i)  is  not  strictly 
correct  at  lower  stresses  (<  1 GPa),  no 
theoretical  analysis  on  the  anisotropic 
contribution  to  wave  profiles  has  been 
made.  Similarly,  there  are  presently  no 
theoretical  or  experimental  studies  as  to 
the  validity  of  assumptions  (ii)-(iv).  The 
good  agreement  presented  here  suggests 
that  these  assumptions  might  be  reasonable 
ones  to  use  for  the  wave  propagation  study. 

Finally,  both  the  experimental  and  the 
theoretical  analyses  presented  here  sug- 
gest different  yield  behavior  for  the 
sintered  and  plasma-sprayed  specimens. 
The  results  on  the  1-D  experiments  and 
the  Hugoniot  data  do  not  agree  with  each 
other  for  sintered  beryllium.  Before 
attributing  this  disagreement  to  the  rate- 
dependent  effect  in  the  Hugoniot  data, 
careful  dynamic  and  static  experiments 
should  be  performed  to  find  out  any  other 
source  of  the  difference.  Future  work 
could  be  concerned  with  the  effects  of 
varying  amounts  of  porosity,  and  with  the 
sintering  conditions  as  they  affect  the 
tradeoff  between  material  strength  and 
its  ability  to  attenuate  a short-duration 
shock  pulse. 
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